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Abstract. We investigate finite volume effects on the pion mass and the pion decay 
constant with renormalization group (RG) methods in the framework of a phenomenological 
model for QCD. An understanding of such effects is important in order to interpret results 
from lattice QCD and extrapolate reliably from finite lattice volumes to infinite volume. 

We consider the quark-meson-model in a finite Euclidean 3-1-1 dimensional volume. In 
order to break chiral symmetry in the finite volume, we introduce a small current quark 
mass. In the corresponding effective potential for the meson fields, the chiral 0(4)-symmetry 
is broken explicitly, and the sigma and pion helds are treated individually. Using the proper¬ 
time renormalization group, we derive renormalization group flow equations in the finite 
volume and solve these equations in the approximation of a constant expectation value. 

We calculate the volume dependence of pion mass and pion decay constant and compare 
our results with recent results from chiral perturbation theory in finite volume. 
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1. INTRODUCTION 

Two important motivations inspire us to study QCD in a finite volume: a generic 
interest in finite volume effects, and the separation of scales afforded by the 
additional length scale. Among nonperturbative methods used in the study of 
QCD, lattice simulations have a very important place. However, they are necessarily 
performed in relatively small volumes, and a good understanding of hnite volume 
effects is required for an extrapolation to inhnite volume. Typically, lattice sizes are 
of the order of L < 2 — 3 fm, while pion masses are of the order m.,^ > 500 — 700 MeV 
[l[ . Thus, we require tools to describe the volume dependence of observables in 

these regimes. In this work, our objective is to describe the volume dependence of 
two of the low energy observables, the pion mass and the pion decay constant, and 
to provide a model for hnite volume effects across a wide range of length scales. 

The low energy behavior of QCD is determined by spontaneous chiral symmetry 
breaking [^: As a consequence, massless Goldstone bosons associated with the 
broken symmetries emerge, and the low energy limit can be described by an effective 
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theory of these light, weakly interacting degrees of freedom. Such descriptions 
in terms of light degrees of freedom only become even better in hnite Euclidean 
volumes: contributions of heavier particles are suppressed by e~^^, where M is the 
typical separation of the hadronic mass scale from the Goldstone masses. Therefore, 
a second motivation to study hnite volume QCD is that it allows isolation of the 
low-energy behavior and a description in terms of the Goldstone modes. 

Different approaches based on this idea have proven very fruitful: Ghiral per¬ 
turbation theory makes predictions for the volume and temperature dependence 
of pion mass, pion decay constant and chiral condensate |^, la, |6| , hnite volume 
partition functions and random matrix theory jl, predict eigenvalue spectra, 
and the dependence of the chiral condensate on volume and quark mass )in| . 

The most important results in the present context are those from chiral perturba¬ 
tion theory (chPT), which relies on an expansion in terms of the three-momentum 
IpI and the pion mass which are small compared to the chiral symmetry break¬ 
ing scale dvr/TT. Gonsequently, a hnite volume places constraints on the expansion, 
since the smallest momentum is determined by the volume, ~ and requires 
diherent expansion schemes, depending on the relative size of L and 1/m^. 


A very useful result obtained by Llischer [11[ relates the leading corrections 


to the pion mass in hnite Euclidean volume to the TTTr-scattering amplitude in 
inhnite volume. The relative shift R[m.,r{L)] of the pion mass mT^{L) in hnite volume 
compared to its value m^(cxo) in inhnite volume is according to this result given by 


R[m^{L)] = 


m,r(-h) —m,r(c>o) 
m,r(c>o) 

3 1 1 

167r2 m^L 


dyF{iy)e-^/^^^^ + 0{e-^^). ( 1 ) 


F{s) is the forward vrTr-scattering amplitude as a function of the energy variable 
s, continued to complex value s, and the sub-leading corrections drop at least 
as where rh > New results have been obtained by using a 

calculation of the vrTr-scattering amplitude in chPT to three loops [nnlo] as input 
for Liischer’s formula jT3.IT^. The shift above the leading one-loop result from this 
approach can then be used to supplement the full one-loop calculation for the 
mass shift in chPT. 

However, this explanation relies purely on the “squeezing” of a “pion cloud” 
whereas - albeit smaller - shifts of meson and hadron masses also appear in 
quenched lattice calculations, show a p ower law behavior for small L rather than 
the behavior expected from chPT |l4i Il5l | , and are generally underestimated . 
This suggests that in addition to the pion effects, quark effects at higher momentum 
scales also contribute directly to the hnite volume mass shifts. 

In the present investigation, we thus use a phenomenological low-energy model 
which is suited to a description of dynamical chiral symmetry breaking and which 
incorporates heavy constituent quarks. As a drawback, this model is not a gauge 
theory and the constituent quarks are not conhned, but decouple from the dynamics 
at low momenta only because of their large masses. Since spontaneous symmetry 








breaking does not occnr in a finite volume, the inclusion of a small quark mass, 
which explicitly breaks the chiral symmetry, is essential. 

We employ an RG method in order to cover the relevant range of length scales 
and pion masses. Thus, we do not have to rely on either the box size or the pion 
mass as an expansion parameter: the RG flow equations remain valid as long as 
the model used as input does. The great advantage is precisely that the RG flow 
equations describe the connection between different momentum scales, and in the 
present case also the dependence on the additional scale 1/L introduced by the 
finite volume. 

2. RENORMALIZATION GROUP FLOW EQUATIONS 

In the present work, we use the two-flavor quark-meson model, an 0(4)-invariant 
linear sigma model where the meson fields (j) = {a, tt) are coupled to two constituent 
quarks through an S'17(2) ® S't/(2)-invariant interaction term. It is an effective 
model for dynamical chiral symmetry breaking below a scale A ~ 1.5 GeV, where a 
description of QGD in terms of hadronic degrees of freedom is valid. Although the 
linear sigma-model by itself is not compatible with the low-energy Trvr-scattering 
data, due to the presence of quarks the low energy constants of chiral perturbation 
theory are reproduced Q 

We consider the model in a four-dimensional Euclidean volume with compact 
Euclidean space and time directions. The bare effective action at the scale A is 

TaM = j d'^x!^q^-dq + gq{mc + (T + if-7i'y5)q+^{df,(j)f + UA{(f))'^, (2) 

where the quark mass term grric explicitly breaks the chiral symmetry. At the scale 
A, the meson potential can be characterized by the values of two couplings: 

Ga(0) = + (3) 

In Gaussian approximation, the one-loop effective action for the scalar fields is 

r[0] = rA[0]-Trlog(rg)[0])+lTrlog(rg)[0]), (4) 

where and r^'^[0] are the inverse two-point functions for the bosonic and 

the fermionic fields, evaluated at the expectation value of the mesonic field, 0. We 
consider an approximation where the field (j) is constant over the entire volume 
and the effective action is reduced to an effective potential, = J d'^x U{(j)). To 
regularize the functional traces, we use a Schwinger proper-time representation 
of the logarithms. The dependence on a cutoff scale k is introduced through an 
infrared cutoff function /a(T/c^) of the form 

k^fairk'^) = -^{Tpy^^exp{-Tk'^) 
ok al 


( 5 ) 


which satishes all required regularization conditions 17, For reasons of con¬ 


vergence of the momentum integrals, we use a = 2. Such a cutoff function makes 
it possible to systematically integrate out only those quantum fluctuations around 
the expectation value which have momenta above the scale k. 

We now wish to obtain a renormalization group flow equation for the effective 
action, which describes how the couplings in the action change with a change of 
the renormalization scale k. We arrive at such an equation by taking the derivative 
of the regularized expression for the effective action with respect to this scale fc, 
and then replacing the bare two-point function from the original expression by the 
renormalized two-point functions, which contain the scale-dependent couplings: 


d 1 dr 


Jo r 


exp[-rr^"y0]] 


k§-J.(re) 


( 2 ) 

B,k 

( 2 ) 


exp[-rr^/;-[0]]. 


( 6 ) 


Note that boundary conditions for the helds are implicit in the traces. In Euclidean 
time direction, the sum over bosonic and fermionic Matsubara frequencies can be 
done analytically before the zero temperature limit is taken. In inhnite volume, 
the momentum integrations for the spatial directions can be performed as well, and 
the flow equation becomes 


d 

fc —17fc((T,7f^L ^ 00) 


3^ 


AN^Nf A 1 \ 

1 i=i * 


( 7 ) 


In hnite volume, the momentum integrations become sums over momentum modes. 
Choosing anti-periodic boundary conditions for the fermionic helds also in the 
spatial directions, we dehne 



The RG how equation for the meson potential in hnite volume then reads 
J ^ TT ( “'2r'i 3/c® V / AN^Nf ^—\ 1 

^ ^ (P+p| + ikr2)5/2 + 


( 8 ) 

.(9) 


The meson masses are the eigenvalues of [17 ^'(ct, P)]*h the second deriva¬ 

tive matrix of the meson potential. In the presence of a hnite current quark mass 
grric, the quark mass = g‘^[{ao + rricY + 2mc{cr — ao) + + “*^ 0 )] 

tains explicitly symmetry breaking terms. Since symmetry breaking terms appear 
only here and are all of this form, we make for the meson potential the ansatz 

4 [^(4-i)] 

Pfc(cr,P) = X] aij{k){a-ao{k)y{a‘^ + n‘^-ao{kfy. (10) 

i =0 j =0 

It depends only on vf^ since the symmetry of the pion subspace remains unbroken. 














TABLE 1. Values for the parameters at the UV- 
scale A which were used in the numerical evaluation. 
The parameters are obtained by fitting the results of 
the RG evolution to a particular pion mass and the 
corresponding pion decay constant in infinite vol¬ 
ume. The physical current quark mass corresponds 
to gruc- 


A 

[MeV] 

muv 

[MeV] 

A(7V 

[MeV] 

U 

[MeV] 

[MeV] 

1500 

779.0 

60 

2.10 

90.38 

100.8 

1500 

747.7 

60 

9.85 

96.91 

200.1 

1500 

698.0 

60 

25.70 

105.30 

300.2 


3. RESULTS OF THE NUMERICAL CALCULATIONS 


We solve the RG flow equations numerically for infinite and finite volume. In order 
to do this, we insert the ansatz m into the flow equations and obtain 

a set of coupled, ordinary differential equations for the couplings aij{k) and the 
minimum (To{k). We choose initial values for the couplings at the scale k = A and 
evolve the couplings until at fc —> 0 all quantum fluctuations have been integrated 
out. Pion mass and pion decay constant are extracted from the resulting effective 
potential. For the case of finite volume, the sums over the momentum modes need 
to be truncated at a finite number of modes. For the results presented here, where 
we use Umax = 40, the effects of this truncation are negligible [21^ . 

In TAB. Q we give an overview over the three parameter sets we have used 
to obtain results for pion masses of 100, 200 and 300 MeV. These parameters are 
obtained by fitting so that the results of the RG evolution match a particular value 
of the pion mass m 7 r(cxo) and the corresponding value of the pion decay constant 
fn{oo) from chPT in infinite volume. With these fixed parameter values, we solve 
the finite volume RG equations, which gives a prediction for the volume dependence 
of m^(L) and f-K^L). The coupling g does not evolve in our approximation and was 
set to g = 3.26, which leads to a reasonable constituent quark mass for physical 
values of A We present our results for the volume dependence of the 

pion mass in FIGS. fTll3l and numerical values for select volumes in TAB. 12 We will 
focus on these results for the pion mass, since they are more easily accessible on 
the lattice than the pion decay constant. For the results on fn{L), we refer to |^ . 

FIGS. fniHl show the relative change i?[m.,r(-h)] = mass 

in finite volume compared to the value mT^{oo) in infinite volume as a 

function of the volume size L. We plot the results for the pion masses mT^{oo) = 
100, 200, 300 MeV on a logarithmic scale. For comparison, we plot in addition 
the one-loop chPT result of j^, and the one-loop chPT result with corrections to 
three-loop order obtained with the Liischer formula 3 . 


The relative change of the pion mass decreases with increasing volume size. 
The RG results are consistently above the results from chPT. Higher loop order 
corrections in chPT increase the mass shift predicted by chPT, and the difference 
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FIGURE 1. Volume dependence of the pion mass for a value ttIt^^oo) = 100 MeV in infinite 
volume. We plot the relative shift of the pion mass from its infinite volume limit = 

(m^(L) — m,r(oo))/m,r(oo) as a function of the size of the volume L. For comparison to the RG 
result, we also plot the results from chPT calculations taken from 0 for the full one-loop result 
(IL chPT) and with nnlo-corrections (IL chPT -|- (nnlo-lo)). 



L[fm] 


FIGURE 2. Volume dependence of the pion mass for a value ttIt^^oo) = 200 MeV in infinite 
volume. For a detailed explanation, see caption of FIG.^ Note the different scales on the axes 
for different values of 


to the RG result becomes smaller. We note that the relative size of the higher 
loop order corrections to the one-loop chPT result become larger with larger pion 
mass. This is not unexpected, since Liischer’s formula becomes an increasingly 
better approximation with increasing pion mass for a given volume size. For 
larger volumes, the differences between the RG result and the one-loop chPT 





FIGURE 3. Volume dependence of the pion mass for a value m,r(oo) = 300 MeV in infinite 
volume. For a detailed explanation, see caption of FIG.^ Note the different scales on the axes 
for different values of m,r(oo). 


TABLE 2. Values for R[mT^{L)], the relative shift of the pion mass in finite volume 
compared to the value in infinite volume (cf. eq. C3l)- We show results for volume 
sizes L — 2.0, 2.5, 3.0 fm, and for pion masses to^(oo) = 10(F 200, 300 MeV. The RG 
results are compared to the exact one-loop chPT results of Q (IL chPT), and to the 
exact one-loop calculation with corrections in three-loop order obtained with chPT 
and Liischer’s formula (IL chPT -|- (nnlo-lo)). In the last column, the difference 
Ai? between the RG result and the three-loop corrected chPT result is given. 


L 

m7r(oo) 

TOttL 


i?[TO,r(L)] 


AR 

[fm] 

[MeV] 


RG 

IL chPT 

{nnlo-lo) 


2.0 

100 

1.01293 

26.6 X 10-2 

8.74 X 10-2 

11.6 X 10-2 

15.0 X 10-2 


200 

2.02586 

5.38 X 10-2 

2.00 X 10-2 

3.31 X 10-2 

2.07 X 10-2 


300 

3.03879 

1.70 X 10-2 

0.56 X 10-2 

1.12 X 10-2 

0.58 X 10-2 

2.5 

100 

1.26616 

10.37 X 10-2 

3.85 X 10-2 

4.97 X 10-2 

5.40 X 10-2 


200 

2.53233 

1.95 X 10-2 

0.73 X 10-2 

1.17 X 10-2 

0.78 X 10-2 


300 

3.79849 

5.31 X 10-3 

1.65 X 10-3 

3.27 X 10-3 

2.04 X 10-3 

3.0 

100 

1.5194 

4.94 X 10-2 

1.91 X 10-2 

2.41 X 10-2 

2.53 X 10-3 


200 

3.03879 

7.85 X 10-3 

2.95 X 10-3 

4.65 X 10-3 

3.20 X 10-3 


300 

4.55819 

1.76 X 10-3 

0.54 X 10-3 

1.05 X 10-3 

0.71 X 10-3 


calculation with the corrections obtained with Liischer’s formula drop exponentially 
as exp(—cmirL) with c > 0, compatible with the error estimate from Liischer’s 
result. As expected, for large volume the slopes of the chPT results and the RG 
results in the logarithmic plot are the same, since the mass shift is then due entirely 
to pion effects and thus should drop as exp(—In contrast to chPT, the RG 
calculation can also be extended to very small volumes, where the chiral expansion 








is considered unreliable, and describe the transition into the regime where chiral 
symmetry is effectively restored. 

In our model, two effects are responsible for the hnite volume mass shift. Chiral 
symmetry is broken dynamically and light Goldstone modes appear when the 
quark helds form a condensate. The presence of the additional scale 1/L affects 
this condensation. In a hnite volume, the lowest possible momentum mode of 
the fermions is a/Stt/L. For small volumes, this acts as an infrared cutoff, which 
effectively “freezes” the quark helds so that fewer modes can contribute to the 
condensate. The mesonic helds, on the other hand, are much less ahected by the 
hnite volume, since the scale /L imposes only a minimum for the smallest non¬ 
zero momentum mode. Fluctuations due to these light mesonic modes lead to a 
further decrease of the condensate, which in turn leads to the observed increase 
in the masses of the light mesons. These pion ehects dominate for larger volumes, 
which explains the convergence of the RG and chPT results for large L. 


4. CONCLUSIONS 

We have presented a study of hnite volume ehects on the pion mass and on 
the pion decay constant. We have applied renormalization group methods to a 
phenomenological model of low-energy QGD with mesons and constituent quarks 
as degrees of freedom. We solved the resulting RG how equations in hnite volume. 
We stress the importance of taking explicit symmetry breaking into account and 
of considering pions and the sigma as individual degrees of freedom. 

We compared our results for the volume dependence of the pion mass to recent 
results obtained in chiral perturbation theory. While they are compatible for small 
pion masses and large volumes, we hnd from the RG generally a larger mass shift 
than the one predicted by chiral perturbation theory. We explain this for small 
volumes by the ehect of the additional IR cutoh L on the quark condensation, while 
for larger volumes the huctuations due to light pions are the dominating ehect. This 
view is supported by the convergence between the RG and chPT results for large 
volumes. Thus, the RG approach in a model including quarks provides a mechanism 
to qualitatively understand the volume dependence. The RG results remain valid 
even to very small volumes, where the chiral expansion becomes unreliable. 

An important point that still requires careful investigation is the inhuence of the 
choice of the boundary conditions for the fermionic helds in the spatial directions. 
As on the lattice, where anti-periodic boundary conditions for quarks lead to a 
larger hnite volume shift than periodic boundary conditions jl^, we observe a 
dependence of our results on this choice. This will be the topic of future work. 
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